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Oxygen isotopic variations in rainfall proxies such as tree rings and cave calcites from South and East Asia 
have been used to reconstruct past monsoon variability, mainly through the amount effect: the observed ^^O 
depletion of rain with increasing amount, manifested as a negative correlation of the monthly amount of 
tropical rain with its 8^^0, both measured at the same station. This relation exhibits a significant spatial 
variability, and at some sites (especially North-East and peninsular India), the rainfall proxies are not 
interpretable by this effect. We show here that relatively higher ^^O-depletion in monsoon rain is not related 
necessarily to its amount, but rather, to large scale organized convection. Presenting 8^^0 analyses of —654 
samples of daily rain collected during summer 2012 across 9 stations in Kerala, southern India, we 
demonstrate that although the cross correlations between the amounts of rainfall in different stations is 
insignificant, the 8^^0 values of rain exhibit highly coherent variations (significant at P = 0.05). 
Significantly more ^^O-depletion in the rain is caused by clouds only during events with a large spatial extent 
of clouds observable over in the south eastern Arabian Sea. 



Stable oxygen isotopic (5^^0) variations of proxies {e.g., cellulose of annual growth rings of trees^'^ and cave 
calcites^'^) of past precipitation are decoded mainly using the significant negative correlation between the 
amount of monthly rain and its 5^^0, known as the 'amount effect'. This effect is observed^"^ in the tropics as 

typically 1.5%o per 100 mm increase in the monthly precipitation. There is, however, a large spatial variation^, 

because of varying cloud condensation processes and their interaction with ambient atmospheric moisture. The 
preferential removal of heavier water isotopologues (H2^^0, and HD^^O) from the cloud vapour during con- 
densation, and of the lighter isotopologue (H2^^0) during rain evaporation during descent, are known to produce, 
respectively, ^^O (also D)-depleted cloud vapour and ^^0(also D)-enriched precipitation. The ^^O-depleted water 
vapour which reaches the sub -cloud layer along with the intense convective down -drafts subsequently moves 
back to the cloud along with the updrafts. This feedback has been suggested to strengthen the amount effect^. 
Though significant amount effect has been observed in the inland stations of India (e.g., GNIP data for Delhi and 
reported data for Jarsuguda^), recent studies over a few Indian sites under the influence of the Indian Summer 
Monsoon (ISM) rainfall such as Mangalore, Kozhikode and Cherrapunji show either a lack of amount effect or 
significant reversed amount effect^°"^^ (i.e., a positive correlation between rain amount and its 5^^0). 

South western India (Kerala) gets the first rains in India during ISM, from the moist air masses originating in 
the southern Indian Ocean and the Arabian Sea and the rainfall is well modulated by the 30-60 day mode of intra 
seasonal oscillations of ISM^^"^^. 5^^0 of ISM rain here has the advantage of minimal alteration by continental 
recycled water (e.g., soil re-evaporation and transpiration from vegetation), and thus could be a good site for 
examining characteristics of the moisture source and the effect of cloud processes on the 5^^0 of monsoon rain. 
Daily rain was collected at 9 stations (Fig. 1) between ~8 to 13°N, —75 to 77.5 °E and 0 to 850 m above sea level 
(above m.s.l) during mid of May to October, 2012. All the stations are located to the west of the Western Ghats 
(which has a maximum height of —2700 m above m.s.l in Kerala). On an average, 73 events were sampled at each 
station, barring events of less than 2 mm/day rain, with a total of 654 samples from all stations. Seasonal total rain 
during 2012 over Kerala was 24% less than its climatological average value (—2040 mm)^^, but individual stations 
recorded both higher and lower values than their respective climatological means. 

In Fig. 2, time series of rainfall (2a), the corresponding rain 5^^0 (2b) are plotted. The daily rainfall over the 
stations varies from 0 to 263 mm whereas the rainfall 5^^0 varies from 0.3 to — 9.7%o from the onset to the 
withdrawal of the monsoon; the amount weighted average 5^^0 of all the stations is —2.8 ± 1.6%o. There exists a 
large spatio temporal variability in the daily rainfall; whereas the 5^^0 values of rain are significantly coherent at all 
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Figure 1 | Sample locations and the amounts & 8^^0 of rain. 

(a) Locations: 1- Thiruvananthapuram (TVM, 53 m asl), 2- Ponmudi 
(PND, 780 m asl), 3- Idukki (IDK, 700 m asl), 4- Ernakulam (EKM, 6 m 
asl), 5- Thrissur (TCR 12 m asl), 6- Palakkad (erstwhile Palghat) (PKD, 
26 m asl), 7- Nilambur (NBR, 35 m asl), 8- Kozhikode (KKD, 30 m asl), 
9- Wayanadu (WYD, 800 m asl), plotted using hcensed CorelDRAW. 

(b) Total rainfall (mm, vertical bars) during Jun-Oct and the 
corresponding seasonal average 5^^0 (filled circles) with la spread, over 
each station, plotted using licensed Sigma Plot. 

the stations (linear cross-correlation coefficients of rain in different 
stations are mostly insignificant, while those of 5^^0 of rain between 
stations are highly significant, as shown in supplementary Table 1). 
Also, the seasonal average of rainfall 5^^0 in all stations agree within 
one standard deviation, whereas the total rainfall values show a large 
variation, with a general increasing trend towards the North 
(Fig. lb), with the exception of NBR, located leeward. The spatial 
variation of the amount of rainfall is induced by the changing oro- 
graphy. The orographic features of Western Ghats play a significant 
role in modulating the ISM rain, since the mountain ranges are 
perpendicular to the prevailing ISM winds^^. This coherent pattern 
of 5^^0 is consistent with the idea that the rain at all the stations is 
derived from a single vapour mass, when clouds have a large spatial 
extent over the region. When we regress the monthly mean 5^^0 
values of each station with the corresponding monthly values of 
rainfall, we find that there is no significant correlation. This means 
that the 'amount effect' according to the prevailing definition does 
not exist in the region. 

Gadgil (2003) noted that OLR values significantly less than 
240 W/m^ and SST greater than 27.5 °C is indicative of large scale 



organised convection^^. Hovmoller diagram i.e., latitude versus time 
plot of Outgoing Longwave Radiation, OLR) averaged over 65- 
77.5°E, (Fig. 2c) indicates the relevance of organised convective activ- 
ity, in determining the temporal variability of rain 5^^0; relatively 
higher ^^O- depletion is seen to correspond mainly with the presence 
of clouds over larger spatial extent (highlighted by arrows). Four such 
events are observed from the first day of the monsoon season until 
the mid- July. During July 15 to August 5, spatial extent of clouds is 
limited, yet continuous rain is recorded in most stations with a peak 
during July 25*. The average amount weighted 5^^0 of rain during 
this period is —1.3 ± 0.3%o, with no corresponding ^^O-depletion 
with the increased rainfall, shows that the rainfall 5^^0 is independ- 
ent of rainfall amount. From August 6* onwards the movement of 
another large cloud cover estabhshed: 5-6 additional, high ^^O- 
depletion events are then recorded. A linear correlation coefficient 
of 0.69 (P < 0.001) is observed between the amount weighted mean 
d''0 with the OLR over the 5°N-12.5°N and 67.5°E-75°E grid box 
(Supplementary Fig. la) till monsoon withdrawal. Thus the pro- 
nounced signatures of organised cloud activity are well captured as 
high ^^O depletion in rain. 

Thus, if 5^^0 of rain is unrelated to the local rain amount and it is 
related to large scale convection, then 5^^0 values of rain should be 
comparable with large scale rain as measured by TRMM (Tropical 
Rainfall Measuring Mission). The difference in the spatial patterns of 
TRMM derived rainfall and OLR during high and low ^^O-depletion 
rain events (i.e. events not associated with large scale convection) are 
shown in Fig. 3. Intense rains (Fig. 3b) from clusters of clouds 
(Fig. 3d) with a large spatial extent over the eastern Arabian Sea, 
including southern peninsular India, are observed during high ^^O- 
depletion rain events. In contrast, low ^^O-depletion rain events 
(Fig. 3a) are characterized by isolated clouds with smaller spatial 
extent (Fig. 3c). The average OLR values over the region of convec- 
tion (5°N-12.5°N and 67.5°E-75°E grid box) during high and low 
^'O-depletion events are 180 ± 19 W/m^ and 235 ± 19 W/m\ 
respectively. The corresponding average TRMM rainfall values for 
the same are 15 ± 7 mm/day and 2 ± 4 mm/day, respectively. Thus, 
the OLR and TRMM rainfall are significantly different during low 
and high ^^O-depletion events. Presence of clouds with large spatial 
extent leads to the higher ^^O-depletion because (i) fresh condensate 
is enriched in ^^O relative to the vapour according to the Rayleigh 
model, and so the increased condensation (higher rain amount) 
reduces ^^O in the rain, (ii) rain drops re-evaporate to a lesser extent 
at the cloud base due to higher ambient humidity and thus less 
enrichment of ^^O compared to events of lesser rain; and (iii) unsat- 
urated downdraft vapour and re- evaporated vapour (both relatively 
depleted in ^^O) during rain events feed the sub -cloud layer 
vapour^'^^'^^, which will be taken into the convective clouds through 
updraft and thereby the rain is further depleted^ in ^^O. 

Strong convection usually occurs in association with strong mid 
tropospheric convergence of moist air and thus vertical velocity too is 
a good surrogate for convective activity. Back-trajectory analysis to 
the mean location (~10.5°N, 76.5°E) of the sampling stations shows 
(Fig. 4b) an upward motion of air masses during events of higher ^^O- 
depletion, while downward or horizontal wind motion is seen during 
low ^^O- depletion rain events. The average trajectory height during 
low ^^O -depletion events is ~2 km, and air parcels do not move 
vertically up when they reach Kerala; in contrast, during high ^^O- 
depletion events, the average trajectory height is below 1.5 km. This 
enables effective pick up of additional moisture from the sub -cloud 
layer and the warm ocean below, before ascending over the rainfall 
site. Between rain 5^^0 and the average vertical velocity (oo, hPa/s) at 
4 pressure levels between 300 hPa and 600 hPa over a grid box of 
75°E-77.5°E, 7.5°N-12.5°N (Supplementary Fig. lb), a linear cor- 
relation coefficient of 0.47 (P < 0.001) is observed up to the time of 
monsoon withdrawal. Recent studies point out the equal contribu- 
tions of both convective and stratiform clouds in the Maximum 
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Figure 2 | Time series of the amount & 6^^0 of rain and convective activity, (a) Time series of measured amount of daily rain (mm), (b) rain 5^^0 (%o) 
relative to V-SMOW (scale on the right), and (c) HovmoUer diagram of average NOAA interpolated OLR over 65°E-77.5°E. Rain events with b^^O 
below — 3.4%o (amount weighted mean of all b^^O values shown as solid line) corresponding to low values of OLR are shown by vertical arrows. Solid 
parallel lines in (c) depict the latitude range of the sampling stations. The vertical Kne indicates date of monsoon withdrawal. For station abbreviations see 
caption to Fig. 1. Plots produced by licensed Sigma plot and GrADS. 
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Figure 3 | Co-occurrence of large scale convection and isotopically depleted rain. Composite of rainfall (TRMM 3B42) and OLR (NOAA) pattern 
during high ^^O depletion (b,d) and low ^^O depletion (a,c) rain events, (e) Shows the days chosen to plot the composite of high ^^O depletion (blue) and 
low ^^O depletion (red) rain events. Maps plotted using GrADS. 
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Figure 4 | Back trajectory analyses, (a) Three-days back trajectory to the 
mean location of sampling stations (~10.50°N, 76.50°E), derived from 
HYSPLIT model, for low ''O depletion events (red lines), high '^O 
depletion events (blue lines) and rain events during October (green lines), 
the month of monsoon withdrawal. Cyan line shows the back trajectory 
ending on 25* July, when low ^^O depletion events of heavy rainfall 
occurred in the absence of organised convection. The filled circles indicate 
the position of the air mass 72 hours prior to rain, and their colours 
indicate the 5^^0 of rain at the final station (plotted at the initial point to 
avoid clustering at the end point); refer to colour code on the right, 
(b) Change in the height of the air mass en-route during high ^^O depletion 
(blue) and low ^^O depletion events (red), and on during 25* July (cyan) to 
the mean cloud base altitude at the sampling station (2 km). Maps plotted 
using licensed MATLAB. 

Cloud Zone (MCZ) in determining the rainfalP°. Large spatial extent 
of stratiform cloud leads^ to mesoscale subsidence of depleted back- 
ground vapour into the sub-cloud layer. ^^O depletion of stratiform 
clouds by this mixing of environmental air could also be an addi- 
tional reason for the observed variability of rain 5^^0. The ^^O 
depleted remnant vapour after the condensation also reaches the 
sub-cloud layer along with the mesoscale subsidence. During the 
convective updraft, this vapour is also taken up into the cloud^'^ from 
the sub-cloud layer and thus the rain becomes more depleted in ^^O. 
Hence, for large scale convective events, along with the other factors, 
recycling of isotopically depleted vapour from stratiform clouds dur- 
ing the MCZ activity also leads to the rain ^^O depletion. 

During ISM, equatorially asymmetric pressure pattern leads to 
cross equatorial flow, also known as the Low Level Jet Stream 
(LLJ), which connects the Mascarene High in the Southern 
Hemisphere and the Monsoon Trough and transports moisture from 
the southern Indian Ocean to south Asia^\ Intra seasonal variability 
of LLJ is also seen associated with the active/break cycles of ISM as 
per the pressure fluctuation between the two hemispheres^ \ Back 
trajectory analysis (see methods) shows that during days of higher 
^^O -depletion, air mass travels longer distances from maritime loca- 
tions from the south of the equator (Fig. 4a). In contrast, during the 3 
days prior to the low ^^O-depletion events, the air mass moves zonally 
from the western Arabian Sea, as well as the inland locations of 
northern Somalia. This underscores the importance of location of 
vapour source and the track the air mass has followed in determining 
the rain 5^^0. The possible reason for the high ^^O depletion in rain 
could also be due to (i) the difference in ^^O content of the marine 
source (e.g., southern Indian Ocean) and (ii) the cloud process along 
the trajectory before reaching the site of observation. 




Longitude (°E) 

Figure 5 | Spatial correlation pattern of the daily average 6^^0 of rain at 
all stations with the TRMM3B42 rainfall over the south Asian region at 
each 0.25 ° X 0.25 ° grid for the period June to mid Oct. Values statistically 
significant at 0.001 level (R > 0.30) are indicated on the right with an 
arrow. Map plotted using licensed MATLAB. 

An anomalous high rain event occurred during 18-28^"^ July with 
low ^^O-depletion all over Kerala. The rain amount was quite high in 
many stations, particularly in the north. Organized convective activ- 
ity was absent (included in the composite plots. Fig. 2c and Fig. 3a). 
Yet, the location of origin of the air mass (Fig. 4a) does not show 
significant difference (taking into account the associated model 
uncertainty) from those during high ^^O- depletion rain events. Air 
mass back trajectory of 25^^ July shows a descending motion up to 
1.5 km (above m.s.l) till the previous 30 hours and subsequently an 
ascending motion (Fig. 4b). Therefore the dominant vapour source is 
likely confined to the atmosphere over the Arabian Sea (with little 
contribution from south of the equator), though the parcel had an 
apparent cross equatorial flow. Thus, along with the absence of large 
scale clouds during this event, this could also be a reason for the 
observed anomaly. 

The daily amount-weighted all-station average 5^^0 of rain and 
grid- wise (0.25° X 0.25°) daily TRMM-derived rainfall shows a stat- 
istically significant (R > 0.30; P > 0.001) negative correlation over 
the south eastern Arabian Sea, including the Kerala coast (Fig. 5). The 
spatial correlation pattern is similar to the OLR and TRMM rainfall 
pattern during the high ^^O depletion events (Fig. 3b & 3d). This 
reinforces our conclusion that this negative relationship is not due to 
'amount effect' (which is local) but rather due to large scale convec- 
tion and the cloud spread. 

After a dry period from mid September to 8^^ October, rain got 
established again. This was the period of monsoon withdrawal and 
the change in wind direction now brought moisture from Bay of 
Bengal (BoB). The rain 5^^0 was much more depleted. The sug- 
gested^'^^'^^ reasons are: (i) the air mass travelled a longer distance 
from the BoB, effectively depleting the rain ^^O over Kerala, due to 
rainfall extraction en route (the continental effect), (ii) Relatively 
lower 5^^0 of surface BoB due to seasonal runoff from the Ganga- 
Brahmaputra river system. 

To summarise, the presence of large scale cloud bands during ISM 
leads to a marked ^^O-depletion in the monsoon rain. This signature 
of convective activity is well captured in rain 5^^0 rather than in the 
local rainfall amount. This is confirmed by an analysis of the spatial 
pattern of clouds using OLR and the TRMM derived rainfall. Intense 
cloud activity with large spatial extent is associated with a strong 
cross equatorial flow and high moist air convergence over Kerala 
during all the high ^^O-depletion rain events. Analyses of proxies 
of convection such as OLR, vertical wind velocity and the air mass 
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back trajectory show high moisture convergence and heavy convec- 
tion in the southern Peninsular India which can effectively take the 
recycled high ^^O depletion moisture from the sub cloud layer. Thus 
paleomonsoon proxies such as teak tree cellulose form Kerala maybe 
interpretable in terms of the degree of large scale convective activity 
during the summer monsoon rather than the amount of monsoon 
rainfall. 

Methods 

Sampling and isotopic analysis. The protocol proposed by International Atomic 
Energy Agency (IAEA) was used for sample collection (http://www-naweb.iaea.org/ 
napc/ih). Using a funnel and a carbuoy, rainfall was collected everyday at 3.30 am 
(UT) and transferred to leak-proof plastic bottles. Isotopic analyses of the samples 
were done with a Thermo Finnigan Delta V plus continuous flow isotope ratio mass 
spectrometer. CO2 equilibration method was adopted for the sample 6^*^0 analysis^^ 
All of the measurements were performed with la external precision of ±0.1%o. The 
precision and accuracy of the measurements were monitored by using an internal 
laboratory water standard, NARM (a large batch of Narmada River water with a 3^^0 
of -4.7%o relative to V-SMOW). We also participated in the Fourth IAEA inter- 
laboratory comparison for oxygen isotopic analysis of water^"^, and values reported by 
us are in good agreement with the IAEA consensus values. 

Monsoon performance in 2012. India Meteorological Department (IMD) defines 
the onset of ISM as the abrupt increase in the rainfall amount in the south-western 
coast of peninsular India with a threshold amount of rainfall (at least 2.5 mm each for 
two consecutive days) recorded in at least 9 out of 14 stations. A certain minimum 
value for OLR and the westerly wind strength in the domain has also been considered 
for the declaration of the onset (http://www.imd.gov.in/section/nhac/dynamic/mok. 
pdf, accessed on 02/05/2013). In 2012, though the monsoon rainfall had started in 2 
stations by 30*^ May itself, rainfall did not establish well in all the stations until 
5* June, declared as the date of onset. Similarly, although the withdrawal of the 
monsoon was declared on 18*^ October, 2012^*^, continuous rain was recorded even 
after 18* October. 

Data sources. Daily rainfall data from TRMM 3B42 with a spatial resolution of 0.25° 
X 0.25° is utilized for the identification of larger scale spatial pattern of rainfall during 
the monsoon season. NOAA interpolated Outgoing Longwave Radiation (OLR) 
data with a spatial resolution of 2.5° latitude X 2.5° longitude and a temporal 
resolution of 24 hrs is used to identify the extent of clouds and their propagation^^ and 
the MERRA (Modern-Era Retrospective Analysis for Research and Applications) 
vertical velocity data between 600 to 300 mb levels are utilized for the analysis of 
convective activity^^. The HYSPLIT (HYbrid Single- Particle Lagrangian Integrated 
Trajectory)^^ model is used for air mass back trajectory analysis for the 72 hrs prior to 
the rain occurrence. 
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